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A technique for maintaining synchronization between optical PPM pulses and a 
receiver clock by means of a delay -tracking loop is described and analyzed. The tracking 
loop is driven by a doubly stochastic Poisson process that contains information about the 
location of the desired slot boundaries. The slot boundaries are subject to slowly varying 
random delays that are ultimately tracked by the loop. The concept of fractional rms 
delay error is introduced to quantify the effects of signal and background induced shot 
noise on the performance of the delay -tracking loop. 


I. Introduction 

Optical pulse-position modulation (PPM) is an accepted 
signaling format for the transmission of information over 
optical channels. With this coding technique, L bits of infor- 
mation are encoded onto one of M = 2 L PPM words by estab- 
lishing a one-to-one correspondence between the possible 
states of L binary digits and the location of an optical pulse 
among M possible slots. If the slot duration is r seconds, then 
each PPM word requires Mr seconds for transmission. For a 
direct-detection system, maximum likelihood decoding con- 
sists of counting the number of photons in each slot and 
selecting the PPM word corresponding to the greatest count. 
To accomplish this, however, requires accurate synchroniza- 
tion between the received PPM slot boundaries and the count- 
ing boundaries defined by the receiver clock. Synchroniza- 
tion for binary PPM optical signals by means of “early-late” 
gates and decision-directed feedback has been considered 
previously in the literature (Refs. 1, 2). Here we consider a 
novel technique for synchronizing M-ary PPM symbols without 
the use of decision-directed feedback, thus effectively remov- 
ing the dependence of the synchronization subsystem on 


decoder performance. The description and analysis of this 
tracking system is the subject of the following sections. 

II. Tracking System Description and 
Analysis 

The slot synchronization technique described here relies 
on the addition of a small “dead time” before and after each 
optical pulse, in order to derive an error signal proportional 
to the delay error that may exist between the received slot 
boundaries and the receiver clock. (Decoder performance is 
not affected if the pulse intensity at the transmitter can be 
increased to maintain constant average pulse count at the 
receiver.) Thus, r = r p + 2r d , where r p is the duration of 
the optical pulse and T d the duration of each “dead-time” 
interval. 

The detector output can be modeled as a doubly stochastic 
Poisson process, represented by a sequence of randomly occur- 
ring unit impulses whose average rate at any given time depends 
on the total signal and background power collected by the 


24 



receiver. Consider the case where the receiver clock is syn- 
chronized with the received optical data (we do not address 
the acquisition problem here). Doppler effects and trans- 
mitter or receiver clock instabilities tend to introduce random 
delay components between the received data and the receiver 
clock which, if uncorrected, may eventually cause loss of 
synchronization. Assuming that the tracking system main- 
tains synchronism between the received data and the receiver 
clock, transmitter clock instabilities and Doppler-induced rate 
changes can be modeled in terms of the delay process Ar = 
Ar cf + At d , whose components are transmitter clock drift 
and Doppler-induced delay, respectively. Typically, At is a 
slowly varying function of time that can be treated as a 
constant over a great many PPM words. 

A block diagram of the delay-tracking loop is shown in 
Fig. 1. The detector output process jc(f;Ar), which contains 
information about the location of the transmitted data- 
pulses, serves as input to the tracking loop. This is a doubly 
stochastic Poisson process with random rate X(f;Ar) = X s (t;Ar) 
+ n b , where X^(r;Ar) is a random rate due to the signal while 
n b is a constant average count rate due to multimode back- 
ground radiation. The receiver clock generates a squarewave 
g(t;AT n ) of nominal period r and noisy delay estimate At = 
lAr + A T cr . The components of are the loop’s estimate of 
the received delay "At, and an independent drift component 
A T cr due to instabilities within the receiver clock. The input 
to the loop filter is the product of the detector output process 
and clock waveform. Defining the delay error as Ar = At - 
this process can be represented as 

y(t;&r e ) = x(t;Ar) g(t;£") (1) 

where the notation implies that the product contains informa- 
tion about the delay error Ar^. Indeed, for Af -ary PPM symbols 
with dead-time, the time averaged mean value given Ar and 
At^, clearly depends on Ar e : 

C T 

F(Ar e ) - — J E {x{t;A7)g(t;&} dt 

l C T 

- Um Jr J M'-At) g(t;£^) dt 

2 n 

= -jgr- Ar e =aAr e ; |At|< 6 (2a) 

5 -min (v|~r d ) ( 2b > 


where 25 is the extent of the linear range,# is the expectation 
operator, n s is the Poisson rate due to signal when a pulse is 
present, and a = 2 yiJMt, Note that J (Ar e ) does not depend 
on the background induced count rate since the time average 
of n b g{t;Ar^) over an integral number of periods is zero. 
Figure 2 illustrates the behavior of y(A r e ) as a function of 
Ar e over the interval | ArJ < r/2, assuming that T d < t/4 . 

It is convenient to decompose y(r;Ar e ) as 

y(t;AT e ) = 7(Ar e ) + A?(r;Ar e ) (3) 

where w(f;Ar e ) A y(f;Ar 6 ) - y(Ar e ) is viewed as an additive 
noise process, containing the conditioning parameter Ar e . 
The power spectral density of this noise process conditioned 
on A r , 5 n (cj;Ar e ), is derived in Appendix A. There it is shown 
that for co = 0 the noise power spectral density does not 
depend on A r e , but takes the value + n b ), where 

n An t /Mt 

is the average Poisson count rate generated by the signal. 
Thus, in a suitably narrow band of frequencies around the 
origin, the power spectral density of the equivalent noise 
process may be approximated as 

(4) 

If A 7 e varies slowly enough to be considered constant over a 
great many PPM symbols, then y(Ar e ) and 5^(co;Ar e ) may be 
approximated by replacing the infinite limits in their defining 
equations by finite values that encompass regions of con- 
stant Ar e . In that case, we conclude that for frequencies 
near the origin, the process y(r;Ar e ) behaves as though it 
were composed of a slowly varying term proportional to the 
delay error; plus an independent additive white noise process 
n(t) characterized by the spectral level (fr + n b ), as defined in 
equation (4). 

Assuming that 1 Ar 1 « 5, the operation of the tracking 
loop can be described by means of the linear loop model 
shown in Fig. 3 (Ref. 3). Using the differential operator 
p A d/di , loop operation can be expressed in terms of the 
equation 

Q 

At = Ar c + A t d - Ar e = —F(p') [aAr € + n] 

(5) 

At A At - Ar 

c — ct cr 

where G is the gain of the voltage-controlled clock and n = 
n(r;Ar f ) as in Eq. (3) (the time index has been deleted for 
the sake of notational simplicity). Solving for A r e yields 
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At = [1 -H{p)] (At c + At d ) -H{p) (6a) 

H(p)AaGF(p)l(p + aGF(p)) (6b) 

The tracking system should be designed to track A t q and A r Q 
accurately, in which case their contributions to A r 6 can be 
ignored. Under this condition the steady-state delay-error is 
due entirely to the equivalent noise process n(t). Assuming 
that \H(Jco )\ 2 takes on significant values only in an interval 
I co | « 2it/Mt , the variance of the delay -error becomes 


Other meaningful fractional rms delay errors may also be 
defined in terms of (oVr^), namely 



o 2 A var (Ar^) 


1 

2m 2 



\H(j u)\ 2 S n (u) do> 
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~T s - <°> 
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(7) 


where B L is the one-sided loop bandwidth defined by 


2 B, 



\H(ju)\ 2 du 


( 8 ) 


Substituting for a and 5^(0) in equation (7) yields 


and using Eq. (2b), 



It should be emphasized that the above equations are meaning- 
ful only if the inequality (pj 5) « 1 is satisfied. Assuming 
this to be the case, Eqs. (11a) and (11c) show that (aj 5) is 
a monotone decreasing function of r d in the interval 0 < r < 
r/4, while maintaining a constant minimum value of 2{pjT p ) 
over the interval r/4 < r rf < r/2. 


B.M 2 t 2 B.t 2 

(n +n ) = - z - z k (ft + n ) 

2 n 2 K * b 2n s s b 


(9) 


From equation (9) it is clear that the fractional rms delay error 
( oJr p ) depends only on B L , n s and n b> and is independent of 
T p for any fixed value of ^ : 



Note that even in the absence of background radiation ( oJr p ) 
does not become arbitrari ly small, but approaches its quan- 
tum-limited value VBJW. 

It is apparent from equations (6b) and (8) that B L depends 
only on the loop filter transfer function F(p ) and the product 
aG. Since G is independent of the other system parameters, 
constant loop bandwidth can be maintained by adjusting G 
so as to keep aG constant. Thus, B L can be held fixed even 
as n s = ar^/2 varies. The behavior of ipjr p ) as a function of 
7f for various loop bandwidths and background rates is shown 
in Fig. 4. 


In typical applications, the loop bandwidth B L is chosen 
large enough to exceed the significant spectral components of 
A r c and A t d , thus allowing the loop to track these quantities 
accurately. The signal-set dimension M and slot-width r are 
generally determined on the basis of communications require- 
ments, while n b depends on the background environment and 
on the particular set of spatial and frequency predetection 
filters employed. In order to achieve the desired decoder error 
probability, some minimum average signal pulse count nr p 
must be maintained, even in the absence of synchronization 
errors. Under the condition (aJB) « 1, random delay errors 
should have little effect on decoder performance if the decoder 
counts photons over the entire synchronous r-second time 
interval. Thus, a value of (aJS) can generally be found to 
guarantee negligible decoder performance degradation while 
simultaneously assuring linear loop operation: let this value 
be denoted by 7. Again using Eqs. (1 la) and (11c), a reason- 
able value for r d can be determined by calculating (a /r ) 
from the given system parameters, and comparing it to 7/2. 
If ( oJt ) < 7/2, then let 






( 12 ) 
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However, if (o It ) > 7/2, then the desired value of ( 0 J 8 ) 
cannot be realized with the initial set of system parameters, 
and hence at least one of the original parameters must be 
altered to accommodate tracking-loop requirements. Generally, 
the tracking loop will perform well, with dead-times substan- 
tially less than 7 / 4 seconds, when operating in high data-rate 
narrow band communication systems where a suitably small 
value of ( °J T p ) can be achieved. 

III. Concluding Remarks 

A technique has been developed for maintaining synchroni- 
zation between a PPM-modulated optical data-stream and a 
receiver clock. This technique relies on the addition of a 
small “dead time” interval before and after each possible 
pulse location to derive an electrical signal proportional to 
the delay error between the received slot boundaries and the 
receiver clock. It was argued that in applications where only 
spectral components near the origin are of interest, and where 
the delay error remains essentially constant over a great many 


PPM symbols, the driving process could be approximated by 
a slowly varying term proportional to the delay error, plus an 
effective additive noise process due to the inherent random- 
ness of the detected optical fields. Based on the above approxi- 
mations, the variance of the delay error due to both signal and 
background radiation was derived, assuming a linear loop that 
effectively tracked out the slowly varying delay error compo- 
nents. We observed that since linear loop operation typically 
implied rms delay errors on the order of a small fraction of the 
linear region, decoder performance should not be signifi- 
cantly affected if the decoder counts photons over the entire 
T-second slot time. It was shown that in applications where 
this is true, a simple algorithm can often be employed to find 
“good” design values for the dead-time interval in terms of 
the given system parameters, and a quantity the designer 
selects to simultaneously ensure linear loop operation and 
acceptable decoder performance. However, a precise deter- 
mination of the dead-time interval (taking into account 
decoder performance degradation and possible non-linear 
effects within the loop) is beyond the scope and interest of 
this article. 
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Fig. 1. Tracking loop block diagram 
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Fig. 2. Time averaged mean value as a function of delay 
error (r d < r/4) 
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Fig. 3. Linear loop model 



Fig. 4. Fractional rms delay error as a function of average 

signal count for various background rates and loop bandwidths 









Appendix A 


The power spectral density of n(t;Ar e ) can be derived 
using the results of Appendix B. Assume that A r € is held 
fixed. Since n(t;Ar € ) = y(r/Ar e ) -y(Ar e ), it follows that 


With P(co) A F {p s (t )} and il(co) A F {p b (t )}, the Fourier 
transform of X(/VAt) g(t,£r ) restricted to the interval 
[- T f T ), becomes 


= S y (wAT e ) - S-ioj'ArJ (A-l) 


where 


where 


(In At \ 2 

S^uAr e ) = 2tt5(w) 

Consider the case where the detector observes signal plus 
multimode background noise fields. Thus, X(/;Ar) = X (f;Ar) + 
n h where X (/;At) is the random rate function associated with 
the signal while n b is a constant average count rate contrib- 
uted by the background fields. The square^wave generated by 
the receiver clock is denoted by g{t;Ar ) as before. For 
M - ary PPM with dead-time the product X(f;Ar) can 

be modeled over the interval [-T, T) as 

\(tAr) g(t;&) = ^ p s + fi.r- jj 

+ E ( iMr + kr - - ) 

i=-K k = 0 \ ' 

(A-2) 

where K = ( T/Mr ) - 1/2, r is the slot duration, and Mt the and 
duration of a PPM word. Elements of the set {£.} are inde- 
pendent random integers taking on values from 0 to (M - 1) 
with equal probability. Given Ar e , p(t) and p b (t) are defined 
as 


F r {X(,;Ar) *<«£>)) - 

( i=-K 

+MP b ( co)7M e /u,AfT / 2 


(A4) 
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(A-3a) 


Recognizing that 


2 £’je' /U,<fi/ ~ fi " )r { 


/=-x: «=-x: 


?»(') = 


n b -,AT e <t<-+AT e 


(A-3b) 


- ; 0 < / < At 6 , — + At 6 < t < r 


= | 7 (w)l 2 +(2A+ 1) 
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and making use of Eq. (A-4), it follows that 


— £ 
2 T 


l/ r r {X^}( 2 




2T 


lim YT f £ Mt.&'J 2 Mt;AT)dt 
r-oc ~ J..T 


= lim 



X(f/Ar) di 


+ \P s (o>)\ (2K + 1 ) 


Af sin ( -r- 
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(A-5) 
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where 


r? s = nrJMr and |g(^z£r^| 2 = 1 


Note that 


lim 
a:— ~ 


ItMI 2 
(2 K + 1) 


oo 



m = ~o» 
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Making use of the fact that 2T = (2 K + 1) Afr, substituting 
the above terms into equation (B-7) and using Eq. (A-l) 
finally yields 


To show this we first observe that 7(co) is a periodic function 
of co with period 2tt/t. Over any interval / that does not 
contain the points c o m - 2irm/r, the denominator 


S n (oj;Ar e ) 


n +n K + 

s b 


l^(co)l 2 

Mr 



M sin (cor /2) 


is bounded away from zero, hence l7(co)l 2 is bounded, and 
l7(co)l 2 


(2**1) 


= 0; co g/ 

m 


Thus, the integral of | ^ (co) | 2 /(2AT +1) vanishes everywhere 
except within arbitrarily small neighborhoods around the 
points co^ as K grows without bound. However, integration 
over any complete cycle yields 


+ |P(co)+AT^(co)| 2 


2tt 


Mr 
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(A-6) 


Thus, the spectral density of the equivalent noise process 
consists of continuous and impulsive terms. The impulses 
occur at integer multiples of the radian frequency 27 r/r; 
however, there is no impulse at the origin since the contribu- 
tion of y was subtracted out. The value of the continuous part 
at the origin is 


1 


2iT(m + lj2)/T 

gM7MI 2 /(2/:+ d = 

27 r(m "l/2)/r 


_L 

M 2 (2 K + 1) 


K M - 1 0 

s sf 

i=-K £-0 


S n (0-,Ar e ) = n s + n b (A-7) 

independent of A r . Therefore, this must be the value of the 
unconditioned spectral density as well: 


27T 

Mr 


S n (0) = {SJO-AtJ} = n s + n b (A-8) 


independent of K . Therefore, this is the value of the integral 
over arbitrarily small neighborhoods around each point oo m as 
K approaches infinity. 

We also make the useful observations that P { 0) = 0, 
|P (0)| 2 = 4 n 2 At 2 , and 


Since for any |A r I < 5, the continuous part of Eq. (A-6) 
begins to deviate significantly from its zero-frequency value 
only as co approaches values on the order of 2tt/Mt, it can be 
argued that in a suitably narrow band of frequencies around 
the origin (namely, \oo\ « 2ir/Mr) the spectral density of the 
equivalent noise process may be approximated by its zero- 
frequency value. 
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Appendix B 






The power spectral density of the product of a doubly 
stochastic Poisson process and a deterministic function is 
derived. Let 


*V) = £ 5('-' e ) (B-l) 

c 

be a sample function of a doubly stochastic Poisson process 
with random rate X(t). Let #(r) be a deterministic function 
and form the product v(r) = g(f) x(f). Define the power 
spectral density of y(t), S (oj), as in Ref. 2 (Chapter 4): 


Thus, 


E tk {ir r ( W )m= 2Z) f 
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X e 
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v m C' 



dt j dt 2 . . . dt k 
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(B-2) 
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(B-5) 


is the Fourier transform of g(r) restricted to the interval 
[- T, T). Using Eq. (B-l) yields 


Y t (oj) = X g(f a )i ,U \- te[-T, D (B-3a) 
c 


i^MI 2 = X X te[-T,D 

C m 

(B-3b) 


First assume that X(t) is known. For a given X(r) and con- 
ditioned on k arrivals in [- T , T), the joint probability density 
of the unordered arrival times t k = (f t 2 , . . . , can be 
expressed as 






re [- T, T) 


(B-4) 


where 


m T A 



X(f) dt 


The first term accounts for those k elements for which SL-m , 
while the second term is due to the remaining (k 2 -k) ele- 
ments of the double sum. Since for a Poisson process with rate 
m T ,E(k) = m T and E{k 2 ~ k ) = , it follows from Eq. (B-5) 

that 


^ £ 


{|r r MI 2 |fc}j = / l£(r)l 2 X(r) 


dt 


+ |F r ^(r)X(0}| 2 


(B-6) 


In general X(r) is itself a sample function of a random pro- 
cess; hence an additional averaging is required over the sample 
space of X(t). Thus, for doubly stochastic Poisson processes, 
the power spectrum is given by the expression 


S y (oj) = lim E 

T -** ® 


lx; 


|g(?)| 2 \(t)dt 


+ \F T {g(t)X(t)}\ : 


(B-7) 
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